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ABSTRACT 

We study the effects of a metallicity variation on the thermal balance and [CII] tine-structure line strengths in interstellar photon dominated 
regions (PDRs). We find that a reduction in the dust-to-gas ratio and the abundance of heavy elements in the gas phase changes the heat 
balance of the gas in PDRs. The surface temperature of PDRs decreases as the metallicity decreases except for high density (n > 10* cm"^) 
clouds exposed to weak (x < 100) FUV fields where vibrational H2-deexcitation heating dominates over photoelectric heating of the gas. 
We incorporate the metallicity dependence in our KOSMA-t PDR model to study the metallicity dependence of [CII]/CO line ratios in low 
metallicity galaxies. We find that the main trend in the variation of the observed CIFCO ratio with metallicity is well reproduced by a single 
spherical clump, and does not necessarily require an ensemble of clumps as in the semi-analytical model presented by Bolatto et al. (1999). 
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1. Introduction 

The bright line emission from photon-dominated regions 
(PDRs) is one of the key tracers of the star formation activity 
throughout the evolution of galaxies in the course of the cosmo- 
logical evolution. Hence, proper modeling of PDR emission is 
of central importance for the interpretation of the observations, 
in order to derive the physical parameters and the chemical 
state of the ISM in external galaxies. The extensive literature 
on PDR emission, both observationally and from the modeling 
side, has largely concentrated on bright Galactic sources and 
starburst galaxies. The effect of different metallicity for the re- 
sulting PDR emission has, up to now, drawn little attention. It 
is, however, very important in order to cover the full course of 
galactic evolution, starting with low metallicity material of cos- 
mological origin. Many nearby galaxies, such as dwarf galax- 
ies, irregul ar galaxies and the Magellanic Clouds have a low 
metallicity (Lisenfeld & Ferrara', 'l998;'Kunth & Osdin, '2000; 
[Pustilnik et al.. .2002: Lee et al., 200 3). Within the Galaxy, as 
well as in other spiral galaxies, there is a radial decrease 
in the metallicity of molecula r clouds and assoc iated HII 
regions (Zaritskv et al-. 1994: | Ari moto., Sof ue & Tsui imotol 
ll996HGiveon et alll2002HBresoUn et al.Ll2004 . These systems 
provide the opportunity to study star formation and photon- 
dominated regions (PDRs) for a variety of metallic! ties. 

In PDRs the molecular gas is heated by the far-ultraviolet 
(FUV) radiation field, either the strong FUV radiation in the 
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vicinity of hot young stars, or weak average FUV fields in 
the Galaxy. The gas cools through the spectral line radiation 
of atomic and molecular species (Hollenbach & Tielens 1999, 
Sternberg 2004). The gas -phase chemistry together with a 
depth dependent FUV intensity lead to the formation of atomic 
and molecular species at different depths through the cloud. 
This typical stratification of PDRs is for example reflected by 
the th e H/H? and C+/C/CO transi tions ( Sternberg & Dalg ariioi 
Il995t iBoger & Sternber i l2005l) . At low visual extinctions 
the gas is cooled by emission of atomic fine-structure lines, 
mainly [CII] 158/im and [01] 63jum. At larger depths, mil- 
limeter, sub-millimeter and far-infrared molecular rotational- 
line cooling (CO, OH, H2O) becomes important together 
with the interaction of dust and gas. Physical conditions 
such as temperature and density can be derived, by com- 
paring the observed line emissions with model predictions 
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[CII] emission is a widely used diagnostic indicator 
of star formation jStacey et all ll99lMPierini et all fl99^ 
iMalhotra et al.'. '2000: 'Boselli. et al.'.'2002':'Pieri ni et alll2003[ 
i Kramer et al., 2004). Observations suggest that low metallic- 
ity systems have higher [CII] to CO rotational line ratios com- 
pared to the Galactic value. In particular, the intensity ratio 
Acii]/^co may vary from ~ 1000 in the inner M iUcy Way, up to 
~ 10^ in extremely low metall i city systems (eg.lMadden et all 
1997: Mochizuki et al., '1998; 'Bolatto et all [1299': 'Madden! 
2000; .Hunter et al...2001.) . Several studies have suggested that 



2 



Rollig et al.: [CII] 158yum Emission and Metallicity in PDRs 



a lower abundance of heavy elements affects the chemical 
stracture of PDRs and the cooling line emission, and that es- 
timates of molecular gas masses from the observed CO(J=1-0) 
line intensities using the standard convers ion facto r may under- 
estimate the true masses in such obje cts j|Wilsoni fr995: Israel, 
11997 '. Israel et al., 2003: Rubio e t allboolk 

IBolatto et al.. (.1999) modelled the metallicity variation of 
the line ratio [CII]/CO(1-0), for an ensemble of spherical 
"clumps", assuming an inverse relation between the size of 
the region and the metallicity. However the sizes of the 
C^, C and CO regions also depend on the chemistry in 
PDR s and the chemical n etwork is modified at low metallic- 
ities jLeaueux et al.lll994l) . Additionally it has been suggested 
that the size of the C^, C and CO regions may also signifi- 
cantly dep end on the o verall cloud morphology, e.g. density 
variations jHegmann & Kegel .2003.) and velocity fluctuations 
llRonigetalll2002h . Moreover the temperature of the molecu- 
lar gas might depend o n the metallicity whi ch affects the ob- 
servable line intensities jWolfire et allll99'5 l). 

We study the effects of metallicity changes on the temper- 
ature and chemical structure of PDRs. In §|2lwe consider the 
dependence of the PDR gas temperature on the metallicity us- 
ing a simplified semi-analytic model and compare it with nu- 
merical results from full PDR model calculations. Our com- 
putations were carried out u sing an upda t ed ve rsion of our 
spherical KOSMA-r model jStorzer et all Il996l) which was 
originally adapted from the plane-parallel model presented by 
Sternberg & Dalgarno (1995). In §|3lwe examine the predicted 
size of the C^ zones as a function of metallicity. We then model 
the strength of the [CII] emission and investigate the the de- 
pendence of the [CII]/CO(J=1-0) line ratio on the metallicity. 
Finally we compare the results with observational data in §|4] 

2. Metallicity dependence of the surface 
temperature 

The basic cooling and heating processes in PDRs, are af- 
fected by the abundances of elements as well as the con- 
tent and the composition of dust grains (I Wolfireetallll995t 
[Kaufman et al- ..1999J . The dust-to-gas ratio {DIG) and the 
optical properties of the dust may depend on the metal- 
licity, Z. Fits to observations suggest that the ratio de- 
pends almost linearly on the meta llicity, DjG oc z' '"^^ 
( iBoselh. Leaueux. & GavazzH l2002l) . There are other stud- 
ies that find deviations from linear ity for hi gher values of 
DIG £isenfeld & FeiTara. L99_8.) . Li & Draine (.2002.) sug- 
gested that the mixture of PAHs in the metal-poor SMC dif- 
fers from the Milky Way. There are a few observations indicat- 
ing that PAHs could have been destroyed by intense UV field s 
at low metalhcities ( iThuan et all Il999t iBolatto et all l2000l) . 
but the detailed composition of dust in low metallicity envi- 
ronments and the influence on its optical properties is not yet 
understood. Because of the insufficient knowledge we assume 
in our model that the composition of the grains does not change 
with metallicity and that the dust-to-gas ratio and the gas-phase 
abundance of heavy elements scale linearly with Z. 

Changes in Z affect the abundances of major coolants as 
well as the electron densities in PDRs. Additionally, a reduc- 



tion in the dust abundance diminishes the UV opacity, the 
photo electric heating rate, and the H2 formation rate. These 
changes affect the temperature and chemistry in the surface lay- 
ers where C^ is most abundant. 

The dependence of the surface gas temperature on 
Z can be estimated considering the balance of cool- 
ing and heating. The dominant cooling processes de- 
pend predominantly on the total hydrogen gas density n. 
[OI]63jUOT, [CII]158jum emissio n, and gas-grain collisions 
are important cooling processes teurke & Hollenbachl Il983l: 
Storzer et al., 1996). Their relative importance in the dif- 
ferent regimes is discussed in Sect. 12.2.11 The dominant 
heating process depends on the far-ultraviolet (FUV; 6- 
13.6 eV) field and the density. Grain photo-electr ic emis- 
sion ( PE1 jd'Hende court & Lege il IL epp & Dal g-arnol 
[l98?, 'Verstrat e et all il99Q>: Bakes & Tielens, 199^, colli- 
sion deexcita tion of FUV pu mped molecular hyd rogen H? 
( Ster nberg & Dalgarnol Il989t [Burton. Hollenbach. & Tielensl 
il990i) and heating from H2 formation play important roles. 
They are discussed in detail in Sect. l2.2.2l By explicitly consid- 
ering the metallicity dependence of each of these cooling and 
heating processes and identifying the dominant processes in the 
different parameter regimes we will show how the energy bal- 
ance in PDRs depends on Z for a quantitative understanding of 
the PDR surface temperature. 

2. 1. The KOSMA-T PDR model 

In our study we use an updated v ersion o f the spherical PDR 
code described in detail by Storz er et alJ ((1996). Briefly, this 
model solves the coupled equations of energy balance (heat- 
ing and cooling), chemical equilibrium, and radiative trans- 
fer in spherical geometry. The PDR-clumps are characterized 
by a) the incident FUV field intensity given in units of the 
mean interstellar radiation field of Draine (1978), b) the clump 
mass, and c) the average density of the clump, for a radial 
power-law density distribution with index y. We incorporate 
the effects of varying metallicity by varying the assumed abun- 
dance of dust grains and heavy elements. The following pa- 
rameters are multiplied by the metallicity factor Z: (a) the to- 
tal effective FUV dust absorption cross section per hydrogen 
nucleus cr; (b) the photo-electric heating rate; (c) the H2 for- 
mation rate; (d) the metal abundance. We consider a range of 
Z from 0.2 to 1. For Z = 1 we use cr - 1.9 x 10"^' cm^, 
C/H=1.4 X 10' ^, and 0/H=3 X 10~^ as st andard values for 
the local ISM (Holl enbach & Tielensl Il99^ . These values are 
slightly lower than re c ent so lar values of 12 -H log(0/H) ^ 
8.7 by | A splyn(^ ^\ ^1,| l|2Q() 4|l. Fo r a detailed discussion see 
iBaumgartner & Mushotzkyl(l2005l) . We do not consider gas de- 
pletion on grains. 

2.2. Semi-analytic approximations 

The results of the full numerical computations can be under- 
stood and anticipated using some simplifying semi-analytical 
approximations. Here, we focus on the surface temperature of 
the PDR at Ay=0, as we are not interested in the shielding prop- 
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erties of low-metal PDRs but in the thermal behavior of the 
low-extinction region dominating the CII emission. We study 
how the thermal properties respond to altered elemental abun- 
dance. We assume the molecular cloud is sufficiently thick such 
that it absorbs all radiation coming from the backside, leading 
to an emission line escape probability fUj = 0) = 1/2 at the 
surface. 

2.2.1. Cooling 

In the simplified model we include the three main cooling pro- 



cm and roil fine cooling 


and gas-grain collisional 


teurke & HoUenbachl 


1 19831 


IStorzeretallll996'). Gas 



cooling is generally dominated by fine structure emission of 
[CII] and [OI]. Gas-grain cooling starts to contribute signifi- 
cantly for high densities. At n > 10^ cm"^ the coupling be- 
tween gas and dust is strong enough so that the efficient cooling 
of the dust by infrared radiation also provides a major coolant 
to the gas. The total cooling rate per unit volume by these ra- 
diative processes is the sum 
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Fig. 1. The solid and dashed line represent the points inn - x- 
parameter space where Aqi = Acn for metallicities of Z = 1 
and Z - 0.2 respectively. 



Atot = Acn -I- Aoi -H Ag_g (1) 

With the abundance of carbon, oxygen, and dust scaling with 
the metallicity, the total cooling rate is also linear in Z. 

Analytic expressions for the three cooling processes are de- 
rived in Appendix A. In Figure we show the relation be- 
tween the fine-structure line cooling contributions as a function 
of n andx for two different values of Z. The region in parame- 
ter space where [CII] cooling dominates over [OI] is shaded in 
gray. For densities above 10"*^ cm""* [CII] cooling is quenched 
and [OI] cooling dominates. At all lower densities [CII] cooling 
dominates. Every point in Figure ^ corresponds to a different 
equilibrium temperature and resulting ratio Acu/Aqi. The de- 
pendence on;^' in Figure^results from the implicit temperature 
dependence of the fine-structure cooling rates. Grain cooling 
dominates only for densities greater than 10'' cm 
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2.2.2. Heating 

The dominant heating process depends on the FUV field in- 
tensity and density. For high intensities grain photo-electric 
heating dominates. T he rate for this process, given by 
iBakes & Tieleiisl ( 1 1994 . is TpE = 10"2^eGo«H Z erg s ' cm^^ 
where e is the photoelectric heating efficiency and Go is 
the UV intensity in u nits of the Habing field. Following 
bakes &Tieiensl ill 994 the photoelectric heating efficiency is 
given by: 



3 X 10-2 



1 + 2x 10-4(Gori/2/nc) 



(2) 



with «(, being the electron density in cm"^ and T the dust tem- 
perature in K. We set Go = 1.71 x 0.5 x to account for the rel- 
ative factor of 1.71 between the Habing and Draine fields, and 
the fact that at the surface of optically thick clouds radiation is 
incident from a solid angle of 2 tt rather than 47r steradians. In 



Fig. 2. The solid and dashed line represent the points inn -X' 
parameter space where rH;,net - TpE for metallicities of Z = 1 
and Z - 0.2 respectively. 



evaluating the efficiency e we use the analytic expression for 
the electron density derived in AppendixiBl 



0.84 X lO^^nZ 1 + A 1 + 14.4 



f{)J5 



cm 



-3 



(3) 



It is common to express the net PE heating rate as Tpg - 
FpE - Aiec, where Aiec is the cooling rate d ue to electron re- 
combination. We adopt the analytical fit from lBakes & TielensI 
(im): 



Arec = 3.49 X 10"'" r 



30 .7,0.944 



n„nZ 



(4) 
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In dense PDRs a second important heating source is the 
colUsional deexcitation of vibrationally excited Hj. The rate 
for this process can be expressed as 



^H* ^ X P nil, AE f ergs 'cm 



(5) 



with «H, is the density of molecular hydrogen , P - 2.9 x 



10- 



s is the pumping rate for a unit FUV field, AE 



23500 K is the characteristic vibrational transition energy, and 
an efficiency factor / accounting for all processes that may re- 
duce the number of de-exciting collisions (see Appendix B). 
The balance equation for the formation and destruction of H2 
is 



(6) 



where D - 2.6 x 10 ^ s ' is the total dissociation rate in a 
X - ^ FUV field, and R = RqZ is the grain surface H2 forma- 
tion rate coefficient (cm-' s"'). Here, it is implicitly assumed 
that dust grains are always covered by enough H atoms, so 
that the recombination rate is only limited by the number of 
H-dust collisions. We use the standard recombination rate Rq 
byflollenbach & Salpeter ( 1971); H oUenbach et al. (197 1): 



^0 = 3 X IQ-^^ fa S T^'^- cm^ s" 



(7) 



where the accommodation coef ficient, fn and the sticking prob - 
ability S are independen t of Z faoUenbach & McKeeL Il979l) . 
ICazaux&Sp aans.'('2004') show that gas-phase formation of H2 
becomes important for Z < 10"^, well below the minimum 
value of Z we consider here. 

The density of atomic and molecular hydrogen is deter- 
mined by Eq.(|6j and can be written as 



1 



1 +20- 



and = n 



a 



1 +2a 



(8) 



respectively, with a - nR/(xD) being the ratio between for- 
mation and destruction rate coefficients. For example, foix - 1 
and n = 10^ cm"-' more than 99% of the gas at the surface is 
atomic. For Z = 1 and a unit Draine field, a = I and nu - wh, , 
for densities n s; 10'' cm""*. As only R depends on Z, it follows 
that a cc Z and the density n at which the H and H2 densities 
are equal at the cloud surface scales as 1 /Z. 
From equations|5]|5J andQit follows that 



Thj = j Ro Z n nu AE f erg s 'cm 



(9) 



We see that Tn* is maximized when iin — n and all of the 
hydrogen is atomic. The efficiency factor / in Equation |5l is 
largest in the limit of high gas density and low FUV intensity, 
where radiative processes become negligible in depopulating 
the excited vibrational levels compared to collisional deexci- 
tation. For a temperature of T - 100 K, Z - \, x - 1' ™d 
n - crcT^ we obtain F™" s; 1.8 x lO"^'* ergcm"-' s"' with 

/ = 2.8 X 10^1 

This is in good agreement with the results from the numeri- 
cal PDR model shown in Fig. lC.ll in Appendixicl The assump- 
tion of a constant formation rate Rq is valid for: x 10^. A 
higher UV field the dust temperature increases leading to a 
rapid reduction of the accommodation coefficient fa- Hence 



in our calculations the maximum H2 heating rate drops for 
X > 10^ as shown in the bottom plot in Fig. lC.ll in Appendix 

o 

H2 not only contributes to the heating, but cools the gas at 
higher temperatures (Sternberg & Dalgarno, 1989). To account 
for the cooling we define the net heating rate V^l - Fh* - 
Ah, . Using the analytic approximations to the molecular level 
structure of H2 derived in Appendixlcjwe obtain 



XP 



■AE 



2 I Acir+Oeir \ 
\ yn I 

9.4 X \(fi^x 



Ah, = 



/ 1.9x10-'' +4.7xlQ-'°;t- 

n nii^ AE y exp(-AE /kT) 
A + D 

X 



(10) 



yn+A+D 
nnii,9A X 10"'^ 



X 



8.6 X 10-^ H- 2.6x10-" 



r exp(-6592K/r) 

X 



rn4-8.6x 10-7 +2.6X 10-" a- 



(11) 



with a collisional rate coefficient y - 5.4 x 10-'^ Vt s-' cm--'. 
As the molecular constants do not depend on the metallicity, 
only the Z-dependence of nn, changes the H2 de-excitation 
heating. 

The relative reduction of the H2 heating at high radiation 
fields is demonstrated in Figure |2l comparing the PE heating 
and the H2 de-excitation heating for the different parameter 
regimes. We see that at any given density Fp^ exceeds F'^' be- 
yond a certain value, but that this limit increases with the gas 
density. 

Although FpE and Fh^ are the two main heating terms it is 
necessary to account for a third process in order to achieve a 
reasonable approximation of the full energy balance. For UV 
fields X < 10^, and densities n < lO'* cm-^, H2 formation 
heating may contribute significantly. Assuming that each for- 
matio n process releases 1 /3 of its binding energy to heat the 
gas ( Sternberg & Dalgarnol Il989l) . the corresponding heating 
rate is: 



2.4 X 10 ^^RntiH erg cm ^s ' 



(12) 



Summing over all three processes we obtain the total heating 
rate 



PE foim 



(13) 



2.2.3. Metallicity dependence 



An inspection of the heating and cooling functions described 
above reveals their metallicity dependence. TableQlsummarizes 
the scaling relations. The radiative cooling functions are linear 
in Z. The photoelectric heating depends on the metallicity via 
e(Z) X Z. The influence of e(Z) can be neglected as long as 
nix ^ 100. Thus, for low UV fields TpE ~ Z. For higher val- 
ues ofx the efficiency accounts for an additional influence due 
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Z=10" 



Fig. 3. Comparison of the KOSMA-t results (open symbols) 
and the semi-analytic values (filled symbols) for the surface 
temperature against the metallicity. The top panel a) is for an 
UV field of ;t' = 10'', the bottom plot b) is for;^' - 10^, and 
the bottom panel c) is for an UV field strength of x = 1 ■ The 
different symbols indicate different surface densities. 



to the electron density which is proportional to Z for high den- 
sities, and independent of Z for very low densities. The metal- 
licity dependence for densities between 10"'... 10^ cm"^ is not 
trivial. Eq.Q shows that the electron density is linear in Z for 
very high n. In the intermediate range this dependence roughly 
shifts from Z" to Z' . This leads to FpE ~ 7} for high values of 
n and YpE ~ Z for very low densities. 

The recombination cooling depends on the metallicity 
through the electron density, resulting in Arec ~ Z and Z'-^ 





low UV field 


high UV field 




;r< 100 


X » 100 


high 


FpE ~ Z 


FpE ~ Z- 


density 


r ~ z 


Fh; ~Z 


(n > 10*^) 


rHT-form 2" 






Arec ~ Z'-5 






A,o, ~ Z 


A,o, ~ Z 


low 


FpE ~ Z 


FpE ~ Z...Z'"' 


density 


Fh; ~Z 


Fh; ~Z 


(n < 10^) 


rHT-form ~ ^ 


rHT-form ^ 




Arec ~ Z 


A,ec ~Z...Z'-5 




A,o, ~ Z 


Atoi ~ Z 



Table 1. Metallicity dependence of the individual heating pro- 
cesses. 



for low and high densities respectively. The Z-dependence in 
the hydrogen heating (Eg. llOland ll \\ comes from the hydro- 
gen density which depends on the metallicity as Z/( 1 +2Z) for 
high densities and low values of;^', and as Z otherwise. 

As a result we show in Fig. |3]the surface temperature of 
model clouds for a variety of different UV field strengths and 
densities computed from the analytic approximation and from 
the full KOSMA-T PDR model. We covered a parameter space 
ranging from n - lO' . . . 10^ cm"^ and x - lO" . . . 10^. It is 
obvious that the metallicity dependence varies strongly over 
the parameter space. We obtain a good agreement for low and 
high UV fields. Even in the intermediate UV and density range, 
where the quantitative accord is weaker, the qualitative depen- 
dence of Tsurf on Z is well reproduced by the semi-analytical 
model. 

The Z-dependence of the temperature can be understood by 
comparing the dominant net rates of heating F/Z and cooling 
A/Z. At high UV fields, where the PE heating dominates, the 
heating is proportional to Z^ at high densities. For a density of 
n = 10^ the PE heating is ~ Z Vz^TT. Due to the high UV 
field n/x < 100 for all given densities thus the electron density 
influences the heating also for small values of n. If the density 
increases the term 14.4r*''^^/n vanishes and FpE ~ "Z} . This is 
reflected in the slopes of the surface temperature in Fig.|31(top). 
For intermediate FUV fields we find a similar behavior with the 
addition that n/;^^ is < 100 for high densities and > 100 for low 
densities, thus the metallicity dependence shifts from Z^ to Z 
with increasing density. This shift can be seen in the middle 
plot in Fig.|3] 

When H2 vibrational de-excitation heating dominates 
(compare Fig.|2ji, the corresponding rate varies as Z/(l + 2Z), 
hence the surface temperature drops as 1/(1 + 2Z). This is 
shown for n - 10'' cm"^ in Fig. |3] (middle) and for n > 
10* cm"^ in the bottom panel of Fig. |3l which gives the temper- 
atures for a FUV field strength of x - 1 ■ For low FUV fields 
and low densities the temperature is proportional to Z due to 
the PE heating as seen in Fig. [S] (bottom). 



6 



RoUig et al.: [CII] IS^/im Emission and Metallicity in PDRs 



Q. 0.1 
Q 




-□- PDR model, a = 
-O- PDR model, a= 1.5 

- ■ -- analytical approximation, a = 

- • -- analytical approximation, a = 1.5 



0.1 



Fig. 4. The width of the layer is plotted against the metal- 
licity. The open circles represent spherical clump of mass M = 
10^ Mo, density no = 10"^ cm^^ FUV field x ^ ^00, and the 
power-law index of the density profile, y = 1.5, whereas the 
open squares represent a spherical clump of mass of 10^ Mo 
with y - 0. The filled circles and squares represent our analyt- 
ical estimates as explained in Sect.|3] 



The offset between the semi-analytical approximation and 
numerical result forx - 10^ is due to a small contribution of 
additional cooling processes in that parameter range. This in- 
creases the total cooling efficiency and hence the temperatures 
in the full numerical calculations are smaller This also holds 
for n - 10^ cm""* and;^' = 1- In that case the cooling is domi- 
nated by CO line cooling which is stronger than [OI] 63/im and 
also by H2O cooling which is comparable to [01] 63/im. Here 
our initial assumptions are clearly underestimating the over- 
all cooling. Even so this does not change the behavior with Z 
which is well reproduced. 

As there is some debate on D/G we tested as an extreme 
example D/G oc instead of linearity. This mainly changes 
the behavior of the heating rates. The dominant surface cool- 
ing processes do not depend on D/G, but only on the elemen- 
tal abundances, while the heating processes are affected by an 
altered D/G. This leads to a decreased heating efficiency for 
Z < 1, hence the surface temperature is significantly lower if 
we assume D/G oc Z^. In the extreme example of x - 10^, 
n = 10^ cm-^ and Z = 0.2 we find r^urf = 240 K, a factor 
of 4 smaller than for D/G oc Z. The cooling inside the cloud 
also depends somewhat on D/G, since the escape probability 
of cooling lines depends on the dust attenuation. 

3. Variation of thie iayer size witli metaiiicity 

At the surface of the PDR the FUV radiation ionizes almost 
all of the carbon atoms. At larger depths the FUV intensity de- 
creases and carbon recombines and is eventually incorporated 
into CO molecules. Thus, a PDR clump can be subdivided into 
a CO core surrounded by an atomic carbon shell and an outer 
C^ envelope. We examine here, the thickness of the C^ enve- 
lope as a function of metallicity Z. We define the C^ envelope 



thickness as the distance from the cloud surface to the depth 
where the abundances of C and C"^ are equal. 

The dominant reaction channels for the formation and de- 
struction of C^ are 

C + v ^ C++e", 
C+ -H e" ^ C -H V 
and 

C+ H- H2 ^ CHJ + V. 

Assuming, that these are the only reactions which influence the 
C^ abundance, the balance equation for the abundance of C^ 
can be written as. 



(14) 



where / is the photoionization rate, ac and kc are the re- 
combination and radiative association rate coefficients. We as- 
sume that the photoionization rate is attenuated exponentially 
as exp(-pAv) where the factor p - 3.02 accounts for the dif- 
ference in the opacity between visual and FUV wavelengths. 
As Ay is determined by dust extinction it scales linearly with 
Z. In the KOSMA-T model we account for an isotropic FUV 
field, and integrate over An ray angles. The ionization rate is 
then given by, 

exp(-3.02Ai/yu) 



3 X 10""V 



Ji 



(15) 



The integral is the second order exponential integral 
£2(3.02 Ay) where fj. - cos0 and is the angle between the 
ray and the normal direction. 

We have defined the radial point rc+ as the location where 
the abundances of C^ and C are equal, nc+ - "c- If we ne- 
glect the contribution of CO at rc* then nc+ = Xcn/2 there. 
Results from the PDR model suggest, that the electron density, 
rieiZ) « 2«c+(Z), thus «c(Z) » nc(Z = 1)Z. For the sake of 
simplicity we chose this expression for rather than the one 
introduced in Eq.|3]which would introduce an additional tem- 
perature dependence. With Xq - 1.4 x 10"^ Z and molecular 
hydrogen dominating the gas density, hh, = n/2, we can re- 
solve Eq.^Jfor the density at rc* 

3xlO-'V£2(3.02Av(rc.)) 
ac Z 1-4 X 10-4 -H 0.5 ;tc 

For a given radial density distribution n(r) and FUV field 
Eq.^]can be numerically solved to obtain rc+, or correspond- 
ingly the width of the C^ layer Dc+ = R-rc+- In Fig.0]we com- 
pare the results from this equation with detailed PDR model 
calculations using the KOSMA-t model for no = 10"* cm'^ and 
X - 100. Here two types of density structures are used: (a) 
n{r) = no (r /R)-t for 0.2R < r < R, n(r) ^ for r > R, and 
n = no 0.2"'' for r < Q.2R, with the total cloud radius R; (b) a 
constant density, n(r) = no, or equivalently 7 = 0. 

The squares representing the constant density model show 
that the C^ layer width depends approximately as Z ' ' on the 
metallicity which closely matches with the as sumption of an 
inverse proportionaUty bv lBolatto et al.l (Il999l) . The circles in 
Fig.|4]represent the model by the y = 1.5. They can be numeri- 
cally fitted with function, Dc+ = (1.195-i-5.758Z)-i47pc.Fig.|4l 
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12+log(0/H) 



Fig. 7. The intensity ratio /[cii]i58/jm/^co(i-0) is plotted against the metallicity (Z = 1 is equivalent to 12+log(0/H)=8.48). The 
lines denote the KOSMA-r results for two different clump masses M - 1 and M - 10 Mq. The density at the surface of the 
clumps and the UV field strength in units of the standard Draine field are given in the plot. The observed ratios of nearby galaxies 
are plotted as squares. 



shows that our results using Eg. II 61 agrees well with the model 
calculations. However the estimated widths are slightly higher 
than the model calculations, which reflects the fact that there 
are more reactions which quantitatively influence the chemistry 
ofC+. 



4. [CII] emission as a function of z 

The [CII] emission as wefl as the [CII]/CO(1-0) fine ratio 
is typically conside red to be a good tracer of star formation 
dStacev et a l.L ll99ll) . The intensity ratio /[cii]i58/jm//co(i-0), ob- 
served in many nearby low metal galaxies, is higher than for 
sourc e s with solar and supe r-solar me talliciti e s ( Madden et aL , 
1997 | iMochizukiet all Il99 8: Bolat to et all 119991: iMaddenI 
2Q0(i lHunteretalll2Q0ll) . Table H summarizes available line 
ratios and metallicities of nearby galaxies. The corresponding 

numbers for Orion are also given as a Galactic re ference. 

Th is dependence has been modelled by iBolatto et alJ 
ill 9991) assuming that the size of the region scales inversely 



with metallicity, and assuming a constant temperature for the 
gas. 

For our spherical PDR model we compute the surface 
brightness as function of clump mass/radius and metallicity. 
The surface brightness is the projected average intensity 



27r I(p)pdp 



71 



(17) 



of the spherical clump, where /(p) is the specific intensity 
along a ray with impact parameter p (iStorzer et al., 1996). The 
line intensities depend on the thermal and chemical structure 
of the cloud. Additionally, for spherical clouds an effective 
area-filling factor of the emissive region has to be considered. 
Particularly lines, which are formed in central regions of the 
cloud are influenced by this area-filling factor. A good exam- 
ple is the surface brightness of '^CO(l-O). For a density of 
n=10^ cm"^ and low metallicities almost the whole cloud is 
devoid of CO due to photodissociation. A higher density or 
metalUcity results in a larger CO core and a higher surface 
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Table 2. Metallicities and observed [CII]/CO(1-0) line ratios of nearby galaxies and Galactic star forming regions. The derived 
values for gas density and FUV strength are also given if available. 



Object 12+log(0/H) [CII]/CO(1-0) References 



NOC 1068 


9.07 


6000 


14,2,12 


NGC 1156 


8.39 


10000 


4 


NGC 1313 


8.25 


1 1740 


14,7 


NGC 1569 


8.08 


37000 


13,4,6 


NGC 4449 


8.3 


14000 


8,6 


NGC 4736 


9 


2500 


\ A 1 '^ in 

14,2,12,10 


NGC 6946 


9 


900 


14,12 


TP 1 n 


8 1 Q 






IZw36 


7.93 


< 3000 


8 


M83 


9.16 


8000 


14,2,5 


M51 


9.23 


3500 


14,9,5 


IC 4622 


8.09 


34000 


3,4 


Orion 


8.75 


6000 


2,11 


LMC 


8.35 


5600 


1 


30 Doradus 


8.43 


69000 


1 


SMC 


8.03 


13000 


1 



References.— (1) Bolatto et al. 1999; (2) Crawford et al. 1985; (3) Heydari-Malayeri et al. 1990; (4) Hunter et al. 2001; (5) Kramer et al. 2005 
(6) Lord et al. 1995; (7) Luhmann et al. 2003; (8) Mochizuki et al. 1998; (9) Nikola et al. 2001; (10) PetitpasrfeWilson 2003; (11) Simon et 
al. 1997 (12) Stacey et al. 1991; (13) Talent 1980; (14) Zaritsky et al. 1994. 




Fig. 5. The surface brightnesses of the [CII] 158jum (dashed- 
dotted) emission lines from the KOSMA-t results (open sym- 
bols) and the approximation from EQ.( ll9t (filled symbols) plot- 
ted against the metallicity for different densities. The cloud 
mass is M=10 Mq, the UV field strength vs. x - 100, the as- 
sumed central temperature is = 35 K, and A - 4. 



brightness. For species which are mainly emitting at the sur- 
face of the cloud this filling effect is negligible. 

In the prior sections we derived approximate expressions 
for the surface temperature of a PDR as well as for the expected 
depth of the envelope. We can use these approximations to 
estimate the total [CII] surface brightness of the PDR. The local 



Fig. 6. The surface brightnesses of the [CII] 158/im (dashed- 
dotted) and i2cO(l-0) 2.6 mm (solid) emission lines plotted 
against the metallicity for two different densities. The cloud 
mass is M = 10 Mq and the UV field strength is = 100. 

emissivity A[cii](r) from Eq. jA.2> can be used to calculate the 
line integrated intensity of the PDR in the optically thin case 

2 ^[CII]('-)_ 



47r 



-dr erg s ' sr ^ 



(18) 



with the total radius R. From Eq. ( I18> follows the mean surface 
brightness of the cloud 



Acii] = ergs cm sr 



(19) 



To calculate /[cujwe assume an exponential temperature pro- 
file T(r) = Tc + TsuiiE2{AAv) with the the second order ex- 
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ponential integral E2, and an arbitrary fitting parameter A. The 
central temperature Tc and the parameter A are different for 
each set of PDR parameters. For demonstration purposes we 
chose = 35 K and A - A Xo estimate the surface bright- 
nesses for multiple PDRs simultaneously. In Fig.|5]we compare 
these approximations with the detailed KOSMA-r results for 
I[cn]- The metallicity dependent behavior is reproduced very 
well and the quantitative agreement is within a factor of 2 as- 
suming the same temperature profile for all 4 models! If we 
drop this assumption and use individual temperature profiles 
for each model the agreement is 10-30%. 

Importantly, the total surface brightness does not scale lin- 
early with the surface density of the clouds. Rather it peaks for 
intermediate values of n, depending on Z. This is mainly a ge- 
ometrical effect, which can be understood by some qualitative 
arguments. If we assume that [CII] is optically thin, we see all 

atoms. In the low density case, where ionized carbon fills 
the whole cloud, the surface brightness is then proportional to 
n (y/A), with the volume of the cloud V and the projected area 
A, hence 7[cii] nR. But R oc « since we kept the cloud 
mass constant, and thus we find /[cii] oz rfil^. For higher densi- 
ties the width of the layer decreases faster than it is compen- 
sated by the growing n. The relative thickness of the layer 
becomes very small for higher densities (i.e. if Dc+ << R). The 
surface brightness then is proportional to nAn R^ Dc* /{n R^), 
hence /[cii] ^ nDc+. We observe a reduced surface brightness 
caused by the geometry of the cloud. This is inverse to the com- 
mon area filling effect for optically thick lines, like e.g. CO (1- 
0), where the projected area of the CO core decreases with de- 
creasing d ensity as dem onstrated in Fig.|5] This was also men- 
tioned bv IStorzer et ali UJ96). This means, that even though 
the local emissivity A[cii] scales linearly with n (see Ea. lA.2> . 
this is not true for the total surface brightness. As a second or- 
der effect we also notice a temperature dependence of the local 
emissivity. The differences between our analytical model and 
the detailed temperature structures from the PDR calculations 
are responsible for most of the deviations shown in Fig.|5l 

We use our PDR model calculations to study the metallicity 
dependence of the [CII]/CO(1-0) line ratio. We adopt a density 
of 10'* cm"^, a y = 1.5, and a UV field , y = 10^ similar to the 
values assumed bv lBolatto et al 1 ill 9991) . Their prediction is an 
average over a clump ensemble as a model for the large scale 
emission from the ISM. In contrast, we start here by investi- 
gating the metallicity dependence for a single, typical clump. 
The discussion below shows that this is already sufficient to 
reproduce the observed trends versus metallicity. The detailed 
investigation of the effects of averaging over a clump ensem- 
ble are left to a subsequent paper Fig.|6lsummarizes the results 
for this typical clump. The dominance of the geometrical effect 
is reflected in the almost constant [CII] surface brightnesses in 
Fig. |6l Hence the line ratio [CII]/CO decreases for increasing 
metallicities and tends to be constant for very high values of Z. 

Figure0shows that the trend in the observed ratios in nor- 
mal galaxies can be represented by a single-clump model with 
M - 1 . . . 10' Mq, shown as solid and dotted lines in Fig.Q 
We also plotted model results for different cloud parameters to 
demonstrate how different observations may be explained by 
different local physical conditions, e.g. the higher observed ra- 



tio for the 30 Doradus region can be explained by a similar 
model, but exposed to an UV field of;^' - lO'* (or alternatively 
by a clump of less mass). This is consistent with derived FUV 
strengths for 30 Doradus jKaufman et all 179991) . The peculiar 
source IZw 36, with a very low [CII]/CO ratio at extremely low 
metallicity can be approximated by a model with a lower FUV 
field of ^ 10, consistent with estimations bv lMochizuki et al] 
(1998). 

Our model results for single clumps reproduce qualita- 
tively the results shown by the semi-analytical clumpy model of 
(Bolatto et al., 1999) for a clump ensemble in reproducing the 
trends versus metallicity. Thus we can confirm Bolatto's find- 
ings whan taking the detailed physical and chemical structure 
of the clumps into account. 

From a practical point of view it is obvious that a clumpy 
ensemble of different clouds should be closer to the true lo- 
cal conditions than a single spherical clump, but we find that 
a clumpy approach is not necessary to explain the observerved 
trend with Z. To model the [CII]/CO line ratio of a particular 
source in detail it may of course be necessary to apply a clumpy 
approach. But to understand the general behavior for different 
metallicities it is sufficient to consider a single, typical clump. 

5. Summary 

We study the effects of metallicity variations on the gas temper- 
ature and [CII] emission line properties of spherical PDRs. We 
find that the surface temperature of PDRs at high UV fields 
varies linearly with metallicity. For low UV fields and high 
densities this metallicity behavior of the surface temperature 
is converse, showing an inverse dependence with metallicity 
due to the dominant H2 heating. We introduce a new two level 
FUV H2 heating and cooling function that properly accounts 
for energy losses via vibrational collisional excitations. 

We examine the dependence of the C^ envelope on metal- 
licity and find that its geometrical depth scales inversely with 
Z. This produces a higher [CII]/CO(J=1-0) line ratio at lower 
metallicities. We used the numerical results from the KOSMA- 
T model to study the dependence of PDR emission lines with 
metallicity. The observed variation of [CII]/CO(J=1-0) with 
metallicity can be explained well by a single-clump model and 
it is not necessary to refere to an average over a clump en- 
semble. We conclude that the [CII]/CO(J=1-0) fine ratios for 
sources with differing metallicities do not provide a strong con- 
straint on the clumpy morphology of a molecular clouds. 
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Appendix A: Cooling functions 

The cooling of the gas is dominated by fine structure line emission of [CII] and [OI]. The line cooling rate can always be written 
as A„/ = nuAuiEuiPijui) erg s"' cm"^ where fi is the escape probabilit y, A„; is the transition probab ility and n„ the number of 
atoms in the upper state u and £„/ is the corresponding transition energy jHollenbach & McKeeifl979i) . Below the critical density 
for the [OI] emission n„ - 8.5 x 10"^ (lOO/T/r)"® cm"^, the main cooling is provided by the 158 fim [CII] line. The general 
cooling rate of a two-level system S can be expressed as: 

Ki = ; —s- erg cm s (A.l) 

l + fexp(E,i/kT)(l + "-^) 

The critical density for collisional de-excitation rta - Aui/jui for the [CII] 158 yum transition is 2.6 x 10^ cm"^. The densities 
we are interested in are in the range of 10^... 10^ cm"^. Inserting the numerical values for [CII] and assuming a relative carbon 
abundance of 1.4 x 10"^ x Z we obtain the cooling rate for the [CII] 158 ;um transition: 

2.02 X 10-24 „z _3 

Acii ; TTTT- erg cm s (A.2) 

l + iexp(92/r)(l + if0) 

The general cooling rate of the [OI]63j(im (^f i -^^ f 2) and [OI]146//m (^Pq Pi) transitions, only accounting for transi- 
tions between neighboring levels, are: 

. A rp nry( noiexp(Eoi/T)gin(n+/3nar,oi) \ _^ _^ 

A12 = Ai2Ei2f)Z I erg cm s (A. 3) 

\gon^ exp(£'oi/r)(n +ySncr,oi)fei n + expiEn/T) g2(n + /3ncr,i2)) I 

Aoi = Aoi£oij6Z ( noigon \ ^^.^^^^-3^-1 ^^ 4-^ 

\gon^ e.xp{EQilT)(n+ I3nc,fi\){g\n + exp{En/T) g2(n + /3ncr,i2)) I 

Inserting the numerical values for [OI] and a relative oxygen abundance of 3 x 10 * x Z in the above equation leads to: 
A63^,,„ = 3.15 X 10""'8.46 x 10"^^ ^Z x 



3 X 10-4 « exp(98 K/D 3 « (« + i j^^^0fL_) 



«2 + exp(98K/r)(n + lj^^^^0^){3n + exp(228K/r)5 (« + i^^M^i^)) 

146//m 



A,46m„ = 1.35 X 10-'M.66x 10"^ tZx 



3 X 10-4 „„2 



«2 + exp(98K/r)(n + lj^^^^0^){3n + exp(228K/r)5 (n + i^^Mpl^)) 



erg cm ^ s ' (A. 5) 



erg cm s ' (A.6) 



This leads to the total [01] cooling rate Aqi = Ags^,,, + Ai46^„. In the high density case gas-grain collisional cooling is also 
contributing to the total cooling of the cloud: 

Ag_g = 3.5 X 10-^"* Vr (r - Tgi-ain) Z . (A.7) 
Tgrain IS the dust temperature as given bv lHoUenbach et al] ( Il99ll) . The total cooling rate is the sum of the above individual rates 
Atot = Acii + Aqi + Ag_g (A.8) 



Appendix B: Tlie electron density at PDR surfaces 

At the surface of a PDR atomic carbon and sulfur are ionized by the impinging, unshielded FUV radiation and atomic hydrogen 
is ionized by cosmic rays (FUV ionization of H is prevented by the Lyman limit). Electrons from dust are negligible due to the 
small number density of dust grains. The relative contribution from these electron donors to the total electron density is shown 
in Fig. lB.ll for the low (left) and high FUV (right) case. The values are computed with the KOSMA-r PDR model. The histogram 
shows the increasing importance of atomic carbon as main electron source with increasing density. Nevertheless the additional 
contributions of atomic hydrogen and sulfur are not negligible. Even at very high UV fields, x ~ 10^, a relevant fraction of the 
electrons is generated in the additional ionization processes. At gas densities ~ lO*"^ cm"^ still 16-20% of the electrons stem 
from the ionization of atomic hydrogen and sulfur This fraction increases rapidly with decreasing density, with the result that at 
a number density of ~ 1000 cm"^ only < 40% of the electrons are due to ionization of atomic carbon. We present an analytic 
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X=1,Z=0.2 X=10^Z=0.2 




density [cm^] density [cm^] 



Fig. B.l. The relative contribution of the various electron donors to the total electron density n^, at the cloud surface as computed in the 
KOSMA-T model. The percentages are given for different values of surface density n (columns) at FUV fields^ = 1 (left) and^ = 10* (right) 
and a fixed metallicity of Z = 0.2. 



approximation to the electron density at PDR surfaces which fits the actual behavior quite accurately and which allows an easy 
interpretation of the metallicity dependence. We assume 



(Xc + Xs)nZ + nH+ 



(B.l) 



Xg and Xs are the relative elem ental abundances of carbon and sulfur ( Xc = 1.4 x 10 Xs = 2.8 x lO' MHollenbach & Tieleiil 
(fl99 9): Federman et al.l(ll993h l. Together with the balance equation for hydrogen: 



^ = Oh nu* tie 



(B.2) 



where ( - 2.3 x 10 s ' is the ionization rate due to cosmic rays dSternberg & Dalgarnol Il995l) . and oh — 3.5 
10~'^(r/300/r) s ' is the recombination rate, we get the following expression for the electron density: 



X 

He — — nZ 
2 



1 + Jl + 



r<"5 /6.05 X 10- 



nZ2 



X 



(B.3) 



with X - Xc + X^. This is in reasonable agreement with the numerical results from the detailed PDR calculations but deviates 
by 5 - 10%. The deviations are due to the different net recombination rates when considering all ionized species. To account for 
this changed rate we have fitted the parameter a in Eq. ( IB.2> to match the electron density from the full PDR model. Assuming 
«H+ ~ Ctin (Irie flH and fitting the constant C to the numerical values we obtain 



X 

— nZ 
2 



1 + Jl + 



r<"5 /7.52X 10- 



nZ2 



X 



(B.4) 



The agreement with the KOSMA-r results is shown in Fig lB.2l Please note that in the high radiation case all C and S at the 
surface are ionized so that the ratio of electrons contributed by them directly reflects their abundance ratio. 



Appendix C: H2 vibrational hieating 

An important heating processes in dense PDRs is collisional deexcitation of FUV-pumped H2 molecules ("Sternberg & Dalgarnol 
Ll995). Here we present a two level approximation for the H2 vibrational heating and cooling valid in the parameter range where 
the process plays a major role (see Sect. 12.2.21 . The approximation reproduces the net heating rate computed by SD95 assuming 
transitions among all 15 vibrational levels in the ground electronic state, but neglecting the rotational structure. 

Vibrational cooling reduces the net heating at large gas temperatures (see Fig. lC.ll bottom). The vibrational cooling is most 
effective at low x for which a large H2 density is maintained. With PDR temperatures of typically less than 2000 K (see Sect. 
13 and the energy gap between the two lowest vibrational levels A/sq,! = 5988 K, we can assume that most of the H2 is always 
in the ground (v = 0) level in this regime. Vibrational cooling is thus basically given by colhsional excitation to v = 1 followed 
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Fig. B.2. A comparison of numerically obtained electron densities (symbols) with the analytical expression (lines) in Eg . IB. 41 at Av=0,x = 10* 
Z = 1 for different surface densities. 



by either radiative decay or photodissociation. Using the molecular constants for the lowest vibrational transition we obtain the 
colhsional cooling rate 



A 



H2 



-A£i,() ri.o exp I —pjT- I n «H, 



(C.l) 



with the spontaneous emission rate coefficient A i,o - 8.6x 10"^ s the coUisional r ate coefficient yi o = 5.4-xl0~ '^ s"' cm 
and the standard photodissociation rate for the v = 1 level of Di - 2.6 x 10"^' s"' jSternberg & Dalgarnolll99"5l) . We found that 
the fit can be improved if AEo.i is increased by 10%. 

In contrast vibrational heating is important when the FUV radiation field provides a significant pumping to higher vibrational 
states. Thus we define a separate equivalent two-level system for the heating. It is characterized by the effective coefficients AiSeff, 
'^eif, Teif, and Deff providing the same heating rate as the full 15 level system 



J •>j 



+ [Aj+xDj]l[yjn] 



= "Hz 



(C.2) 



(C.3) 



1 + [Aeif +xDs«]l[y^«n\ 

The quantity P, denotes the formation rate of vibrationally excited H2 for the different levels, Ptot represents the sum rate over 
all levels. The effective coefficients can be easily obtained by considering different asymptotic values of the density n and the 
radiation field This yields Aot ■ AEeff = 9.4 x 10"^^ erg s"', yeff = ri,o,£*eff = 4.7 x 10"'"s"', and Aeif = 1.9 x lO^^s"'. A 
comparison of the heating rates using our effective two-level system with the results using the 15-level molecule in our KOSMA- 
T model is shown in Figure ICTI Here, we scan the parameter ranges where the H2 vibrational heating and cooling gives a major 
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Fig. C.l. The net H2 heating rate is plotted for different values of surface density n at a x = 1 FUV field (top) and for different UV field 
strengths for a fixed density of n = 10^ cm"' (bottom) over the temperature. The symbols are the numerical results from the 15-level system in 
the KOSMA-T model, the curves show the results using our effective two-level approximation. For^ = 1 we also plot a low Z case (solid, Z = 1, 
and dashed, Z = 0.2). The deviations at the highest temperatures and radiation fields are due to cooling transitions from higher vibrational levels 
in this regime which are ignored in our approximation. 
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contribution to the overall energy balance (see Fig.|2ji. In the upper plot representing low radiation fields and varying densities we 
find an almost perfect agreement of the two level approximation with the full numeric treatment. At high densities and varying 
radiation fields shown in the lower plot, we find a good match at temperatures below about 3000 K. The deviation at higher 
temperatures is due to the neglect of cooling contributions from higher vibrational levels. It has no impact on the overall PDR 
model, because photoelectric heating clearly supersedes the vibrational contribution at these conditions. 

With the simple analytic two-level approximation, we can easily understand the quantitative behavior of the H2 vi- 
brational energy bala nce from basic principles thus providing a handy tool for estimates of temperature structures. 
burt on. Hollenbach, & TielensI (Il990) also introduced a two-level approximation for the H2 heating. However, they considered 
only a single pseudo excited level with an energy corresponding to v = 6 and hence did not properly account for cooling via rapid 
excitation to v = 1, or heating via pumping and collisional deexcitation from all 15 levels. 



